Sequence comparison of the mitochondrial genomes of Plesionika species (Caridea: Pandalidae), gene rearrangement and phylogenetic relationships of Caridea

Background Despite the Caridean shrimps’ vast species richness and ecological diversity, controversies persist in their molecular classification. Within Caridea, the Pandalidae family exemplifies significant taxonomic diversity. As of June 25, 2023, GenBank hosts only nine complete mitochondrial genomes (mitogenomes) for this family. The Plesionika genus within Pandalidae is recognized as polyphyletic. To improve our understanding of the mitogenome evolution and phylogenetic relationships of Caridea, this study introduces three novel mitogenome sequences from the Plesionika genus: P. ortmanni, P. izumiae and P. lophotes. Methods The complete mitochondrial genomes of three Plesionika species were sequenced utilizing Illumina’s next-generation sequencing (NGS) technology. After assembling and annotating the mitogenomes, we conducted structural analyses to examine circular maps, sequence structure characteristics, base composition, amino acid content, and synonymous codon usage frequency. Additionally, phylogenetic analysis was performed by integrating existing mitogenome sequences of true shrimp available in GenBank. Results The complete mitogenomes of the three Plesionika species encompass 37 canonical genes, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). The lengths of these mitogenomes are as follows: 15,908 bp for P. ortmanni, 16,074 bp for P. izumiae and 15,933 bp for P. lophotes. Our analyses extended to their genomic features and structural functions, detailing base composition, gene arrangement, and codon usage. Additionally, we performed selection pressure analysis on the PCGs of all Pandalidae species available in Genbank, indicating evolutionary purification selection acted on the PCGs across Pandalidae species. Compared with the ancestral Caridea, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species—P. izumiae and P. lophotes. We constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes. The results revealed that family Pandalidae exhibited robust monophyly, while genus Plesionika appeared to be a polyphyletic group. Conclusions Gene rearrangements within the Pandalidae family were observed for the first time. Furthermore, a significant correlation was discovered between phylogenetics of the Caridea clade and arrangement of mitochondrial genes. Our findings offer a detailed exploration of Plesionika mitogenomes, laying a crucial groundwork for subsequent investigations into genetic diversity, phylogenetic evolution, and selective breeding within this genus.

In the present study, we sequenced and analyzed three complete mitogenomes of Plesionika species (i.e., P. ortmanni, P. izumiae, and P. lophotes).Our objectives were

Sampling, identifcation and DNA extraction
Three wild species of P. ortmanni, P. izumiae and P. lophotes were collected from two different sea areas in Zhejiang Province, China (Table 1).Experts from the Marine Biology Museum of Zhejiang Ocean University morphologically identified the specimens, referencing the literature (Kim et al., 2012;Li, 2006).The three Plesionika species share characteristics such as having no dorsal ridges or protrusions on their abdomens.Additionally, the bristles on the antennal stalk are sharp and pointed, extending to the distal edge of the first antennal segment.The differences include: in P. ortmanni, the sixth abdominal segment's length is 1.5 times its maximum height; the tail fan's length is also 1.5 times that of the sixth abdominal segment; three pairs of spines are present on both the dorsal and posterior margins.The length of the antennal scale is 4.3 to 4.4 times its width.The second pair of walking legs is approximately equal in size.In P. izumiae, the sixth abdominal segment is 1.7 times its maximum height; the tail fan is 1.4 times the length of the sixth abdominal segment; there are three pairs of small spines on the dorsal margin; and there are three pairs of spines on the posterior margin.The length of the antennal scale is 4.2 times its width.The second pair of walking legs is unequal in size.In P. lophotes, the sixth abdominal segment is 1.5 times its maximum height; the tail fan is 1.6 times the length of the sixth abdominal segment; there are four pairs of small spines on the dorsal margin; and there are three pairs of spines on the posterior margin.The length of the antennal scale is approximately 3.4 times its width.The second pair of walking legs is unequal in size.Samples were preserved in absolute ethanol before DNA extraction.The total DNA was extracted using the salt-extraction procedure and stored at −20 • C for sequencing (Aljanabi & Martinez, 1997).

Mitogenomes sequencing, assembly, and annotation
The complete mitogenomes of three Plesionika species were sequenced by using nextgeneration sequencing (NGS) on the Illumina Hiseq X Ten platform by Origin gene Bio-pharm Technology Co. Ltd. (Shanghai, China).The mitochondrial genomic DNA of the samples underwent initial quality control, where 1.0% agarose gel electrophoresis was utilized to assess the quality of the DNA.Additionally, a nucleic acid quantifier (NanoDrop, Wilmington, DE, USA) was employed to detect the purity and concentration of the DNA.The quality-controlled mitochondrial genomic DNA of the samples was randomly fragmented into 300-500 bp segments using a Covaris M220 ultrasonic disruptor.The fragmented DNA was subsequently purified to construct sequencing libraries.The steps involved are as follows: DNA end repair, 3 adenylation, sequencing adapter ligation, and recovery of target fragments through agarose gel electrophoresis.The final sequencing libraries were generated through PCR amplification and were subsequently subjected to sequencing on the Illumina HiSeq TM platform.Trimmomatic v0.39 software was used to filter out low-quality reads, duplicate reads, sequences with a high ''N'' ratio, and sequencing adapter sequences (Bolger, Marc & Bjoern, 2014).The reads from the three species were reassembled using the NOVOPlasty assembly software (Dierckxsens, Mardulyn & Smits, 2016).The assembled sequences were compared with other Plesionika species genomes in GenBank, and the COI and 16S rRNA sequences were verified using NCBI BLAST (Altschul et al., 1997).Aberrant start and stop codons were identified by comparing with similar codons in other invertebrate species.The online software MITOS (Bernt et al., 2013) was utilized for structural and functional annotation, with manual corrections performed to obtain the final complete mitogenome.The sequenced mitogenomes were uploaded to the GenBank database (Table 1).

Sequence analysis
The circular visualization of the mitogenomes of the three Plesionika species was completed using the CGView server (Grant & Stothard, 2008).The analysis of the nucleotide composition of the whole mitogenome, protein-coding genes (PCGs), rRNA, tRNA genes, and AT content was conducted using MEGA-X (Kumar et al., 2018).The base skew values were calculated using the formulas AT-skew = (A − T)/(A + T) and GC-skew = (G − C)/(G + C) (Alexandre, Nelly & Jean, 2005).The accuracy of transfer RNA genes and their secondary structures was confirmed using MITOS (Bernt et al., 2013).Base composition, nucleotide composition and relative synonymous codon usage (RSCU) of each protein-coding gene were calculated using MEGA-X (Kumar et al., 2018).The Ka/Ks ratios of the three mitogenomes were estimated using DnaSP 6.0 (Rozas et al., 2017).

Gene order analysis
For comparative analyses, an additional 79 complete mitogenomes of Caridea were obtained from GenBank (Table S1), in addition to the three mitogenomes sequenced in this study.The gene arrangements of all 82 mitogenomes were compared with the ancestral Decapoda, with the aim of identifying potential novel gene orders that have not been reported in previous studies.To ensure that observed gene order differences were not caused by mis-annotations, any mitogenome in Caridea that deviated from the ancestral pattern underwent re-annotation using MITOS (Bernt et al., 2013).

Phylogenetic analysis
To explore the phylogenetic relationships of Pandalidae, sequences of 79 species from 11 families within Caridea were downloaded from GenBank (Table S1).The mitogenomes of Hemigrapsus sinensis (NC_065995) and Helicana japonica (NC_065158) from Brachyura served as outgroup, and phylogenetic analyses were performed based on the 13 PCGs of these 84 species.The sequences of 13 PCGs from each sample were identified using DAMBE 7 software (Xia, 2018).The PCG sequences of these 81 species were aligned with MEGA-X's ClustalW (Kumar et al., 2018).Subsequently, Gblocks v.0.91bREF was applied to determine and select conservative regions, removing divergent and ambiguously aligned blocks (Castresana, 2000).DAMBE 7 was used to assess the suitability of these sequences for phylogenetic tree construction (Xia, 2018).
The analysis of phylogenetic relationships was conducted employing maximum likelihood (ML) method through the program IQ-tree 2.1.3(Minh et al., 2020) and Bayesian inference (BI) method through the program MrBayes 3.2.7a(Ronquist et al., 2012) to analyze the phylogenetic relationships.The ML tree was constructed using the program IQ-TREE (Minh et al., 2020), with the best substitution model (TIM2 + F + R7) being filtered based on Bayesian information criterion (BIC) via ModelFinder (Kalyaanamoorthy et al., 2017) calculations, and 1,000 ultra-fast bootstraps were set for the rebuilding of the consensus tree.For BI tree construction using MrBayes v3.2, the initial step involved using PAUP 4 (Swofford, 1993) for format conversion.This was followed by the integration of PAUP 4, ModelTest 3.7 (Darriba et al., 2020) and MRModelTest 2.3 (Nylander, 2004) within MrMTgui to identify the optimal alternative model (GTR + I + G) as dictated by the Akaike information criterion (AIC).Four Markov Chain Monte Carlo (MCMC) chains were simultaneously run for 2 million generations, with a sampling frequency of every 1,000 generations.During the initial burn-in phase, 25% of trees were discarded, and convergence of independent runs was evaluated by the mean standard deviation of the splitting frequency (<0.01).Finally, the phylogenetic tree was edited using the software FigTree v1.4.3 (Rambaut, 2018).
Both ML and BI methods were used to construct phylogenetic trees.The non-parametric bootstrap support values and Bayesian posterior probabilities generated by these two methods represent the support rates of nodes, respectively.The non-parametric bootstrap method tends to underestimate support rates of nodes, whereas Bayesian method tends to overestimate them (Suzuki, Glazko & Nei, 2002).A maximum likelihood value greater than 70% indicates that the clade relationship is well resolved.Conversely, a value between 50% and 70% is considered weak support, and anything below is regarded as unresolved (Huelsenbeck & Hillis, 1993).Similarly, a Bayesian posterior probability of 95% or higher indicates that the clade support rate is well established (Leaché & Reeder, 2002).
The nucleotide compositions of these three newly sequenced mitogenomes were: A: 32.70% to 35.91%, T: 31.46% to 31.94%,G: 11.60% to 13.85%, C: 20.55% to 21.98% (Fig. 2A).The contents of A and T exhibited high values, indicating that codon usage was biased towards A and T, which is consistent with the reported complete Pandalidae mitogenomes (Sun et al., 2020).The three species had low G and C contents, indicating obvious bias against G and C.Among them, the AT contents ranged from 64.16% to 67.85%, while the AT-skew were positive in the range of 0.019 ∼0.058 and GC-skew were negative in the range of −0.227 ∼−0.285 (Fig. 2B).

Protein-coding genes and codon usage
The PCGs in these three Plesionika species mitogenomes had total length of 11,192 bp (P.ortmanni), 11,134 bp (P.izumiae) and 11,041 bp (P.lophotes), respectively, including seven NADH dehydrogenases (ND1-6 and ND4L), three cytochrome oxidases (COI -III ), two ATPases (ATP6 and ATP8) and one cytochrome b (Cytb) (Fig. 1, Table S2).Among these species, the ND5 was identified as the longest PCG, ranging from 1,644 to 1,719 bp, while the ATP8 gene consistently presented as the shortest, with a uniform length of 159 bp.High AT contents were also observed in the base composition of these species, with P. lophotes exhibiting the highest AT content at 66.05%.Additionally, the AT-skew values were found to be negative, ranging from −0.170 to −0.178 (Fig. 2).Upon comparison of initiation and termination codons of all PCGs of the three Plesionika species, we found five initiation codons and two termination codons.The PCGs of these three mitogenomes were predominantly initiated with ATG, ATT, and ATA.Exceptions included the ATP8 of P. ortmanni and ND3 of P. lophotes, which started with ATC, and the COI of P. ortmanni, which commenced with ACG (Table S3).The majority of PCGs in these three mitogenomes terminated with TAA and TAG, while ND4 in all three mitogenomes, along with COI in P. ortmanni and ND5 in P. lophotes, terminated with a single T. The occurrence of incomplete termination codons is a notably common phenomenon within the mitochondrial genes of both vertebrates and invertebrates (Hamasaki et al., 2017).
The analysis of amino acid compositions in PCGs in the three newly sequenced Plesionika species were relatively similar (Fig. 3, Table S4).The most frequently used amino acid is Asn, followed by Leu1, Lys, Phe, Pro, and Thr, while Arg and Cys are less commonly used amino acids.Comparing the relative synonymous codon usage (RSCU) of 13 PCGs in the three species, the result showed that the usage frequency of UUA (Leu), UCU (Ser) and AUA (Met) codons in their sequenced mitogenomes was higher.In P. ortmanni, the highest RSCU was found for CGA (Arg), followed by UUA (Leu), ACU (Thr), UCU (Ser) and AUA (Met).In P. izumiae, the highest RSCU was found for UUA (Leu), followed by AUA (Met), CAA (Gln), UCU (Ser) and CCU (Pro).In P. lophotes, the highest RSCU was found for CGA (Arg), followed by UUA (Leu), UCU (Ser), CCU (Pro) and AUA (Met).The lowest RSCU in all three species was observed for GAG (Glu).

Transfer and ribosomal RNAs
In common with other Caridea mitogenomes, the mitogenome of three Plesionika species contains 22 tRNA genes (Fig. 1, Table S2).The total length oftRNAs in three Plesionika species mitogenomes were 1,471 bp (P.ortmanni), 1,465 bp (P.izumiae) and 1,476 bp (P.lophotes), and the length of tRNAs in these species ranging from 59 to 72 bp (Table S2).All tRNAs exhibited high AT contents, with the AT content for the three species being 66.96% (P.ortmanni), 66.21% (P.izumiae), and 67.48% (P.lophotes) (Fig. 2A).The tRNA genes of P. ortmanni had a weakly negative AT skew (−0.005) and positive GC skew (0.115), while the tRNA genes of P. izumiae and P. lophotes had positive AT skew (0.021 and 0.034, respectively) and GC skew (0.083 and 0.079, respectively) (Fig. 2B).The examination of the secondary cloverleaf structure of the 22 tRNAs from these species was conducted.In P. ortmanni, it was found that trnS1 could not form a secondary structure due to the absence of dihydrouracil (DHU) arms, a phenomenon commonly observed in metazoans (Yamauchi, Miya & Nishida, 2003) (Fig. 4).Additionally, it was noted that the trnP gene in P. ortmanni, the trnA gene in P. izumiae, and the trnD, trnF, and trnH genes in P. lophotes were lacking the T C loop.In contrast, the remaining genes displayed a typical cloverleaf structure (Rich & Rajbhandary, 1976).Comparing the tRNA genes of the three species, it was found that each corresponding amino acid was encoded by the same anticodon.The total lengths of the 16S rRNA and 12S rRNA genes were similar in three species, with P. ortmanni, P. izumiae, and P. lophotes having total lengths of 1,361 bp, 1,328 bp, and 1,322 bp for 16S rRNA, and 803 bp, 811 bp, and 812 bp for 12S rRNA, respectively (Table S2).The 16S rRNA and 12S rRNA genes of three species were located between trnL1 and trnI, and were separated by trnV.High AT contents were also demonstrated, with the AT content for the three species being 67.83% (P.ortmanni), 69.98% (P.izumiae), and 70.76% (P.lophotes) (Fig. 2A).The rRNA genes of P. ortmanni had a weakly positive AT skew (0.012) and positive GC skew (0.253), while both the rRNA genes of P. izumiae and P. lophotes had negative AT skew (−0.059 and −0.056, respectively) and positive GC skew (0.349 and 0.339, respectively) (Fig. 2B).

Selective pressure analysis
In genetics, the Ka/Ks ratio, which represents the ratio between the nonsynonymous substitution sites (Ka) and the synonymous substitution sites (Ks), is commonly used to understand the dynamic evolution of PCGs.In this study, the Ka/Ks ratios of the 13 PCGs were calculated using the nine sequenced Pandalidae species (Table S1) to investigate the relationship between evolution and selection pressure (Fig. 5).Results showed that the Ka/Ks ratios of PCGs range from 0.088 (COI )to 0.341 (ATP8).The Ka/Ks ratio of the COI gene was the lowest, indicating that the COI gene was under the greatest selection pressure and the gene sequence was relatively conservative.As a result, it is widely used as a potential molecular marker in species identification and phylogenetic studies (Astrin et al., 2016).
In general, a gene is considered to be positively selected when the Ka/Ks is greater than 1, neutral evolutionary when the Ka/Ks is equal to 1, and purified selected when the Ka/Ks is less than 1 (Nei & Kumar, 2000;Yang, 2006).In this study, the Ka/Ks ratios of the 13  S1.Full-size DOI: 10.7717/peerj.17314/fig-5 PCGs genes were less than 1, indicating that the genes of Pandalidae species were subjected to purification selection during evolution.

Gene rearrangement
Mitochondrial gene arrangement is an important tool for the study of systematic geography and phylogeny, which provides informative insights into the evolution among metazoans (Beagley, Okimoto & Wolstenholme, 1998;Searle, 2000).In general, mitochondrial gene arrangement is relatively stable in vertebrates, such as fish, amphibians, and most mammals (Fu, Chen & Zhao, 2009).However, in invertebrates, varying degrees of gene rearrangement are commonly observed in mitogenomes (Ye et al., 2021;Boore, 1999).
The gene orders within the infraorder Caridea mitogenomes were compared with those of ancestral Decapoda in this study.It was found that the mitochondrial gene orders (MGOs) of the families Atyidae, Alvinocarididae, Acanthephyridae (Spence Bate, 1888), Oplophoridae (Dana, 1852), and Nematocarcinidae (Smith, 1884) matched those of the ancestral Decapoda.Meanwhile, gene rearrangement was identified in 24 species across 6 families of Caridea (Fig. 6).This contradicts the previous view that the gene order in the Caridea is conservative (Wang et al., 2018;Miller et al., 2005;Ivey & Santos, 2007;Lü et al., 2019).

Phylogenetic relationships
In this study, we constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes.The analysis included 82 caridean species, focusing on P. ortmanni, P. izumiae, and P. lophotes, with H. sinensi and H. japonica serving as outgroups.The topological structures of the phylogenetic tree reconstructed using two methods are identical (Fig. 7), but there were slight differences in the support values of some of the clade branches.The support values of BI were generally higher than ML, with the majority of the nodes having a support value of 1.On the other hand, the support values of ML, except for one node with a support value of 53 (the node between Alpheus japonicus and Alpheus randalli), were between 78 to 100 for all other nodes.The phylogenetic tree analysis demonstrated that all families within Caridea exhibit monophyly.Strong monophyly was exhibited by the families Pandalidae, Thoridae, Lysmatidae, and Hippolytidae, which clustered together forming a single clade.Acanthephyridae and Oplophoridae were closely related in phylogenetic relationship, forming a sister group, and subsequently clustered with Alvinocarididae.A large clade was then formed by these three families with Nematocarcinidae.Additionally, the families of Alpheidae and Palaemonidae were also found to be closely related, forming a sister group relationship.
The phylogenetic tree revealed that these three newly sequenced species within the family Pandalidae did not cluster together.Phylogenetic analysis showed that P. izumiae and P. lophotes were closely related and formed a distinct branch, while P. ortmanni, P. sindoi and P. edwardsii clustered together on another separate branch.At the genus level, we observed that the genera Plesionika and Heterocarpus within the family Pandalidae were not monophyletic.Two distinct clades were formed by the five species of the genus Plesionika, while Parapandalus sp. was grouped with the genus Heterocarpus.

DISCUSSION
It was revealed by a comparison of the base content of the mitochondrial genomes of three Plesionika species that a higher AT base content than CG was exhibited by all three genomes, a phenomenon commonly observed in the mitochondrial genome sequences of decapods (Sun et al., 2020;Zhu et al., 2021;Wang et al., 2019).The trnS1 gene of the P. ortmanni lacked the DHU arm, preventing it from forming a typical cloverleaf structure.The loss of the DHU arm in trnS gene has been reported in most mitochondrial genome studies of Caridean shrimp (Ye et al., 2021;Shen et al., 2012).Selective pressure analysis was conducted separately on the PCGs of nine species in the Pandalidae family, and the results showed that the Ka/Ks values of all 13 PCGs were less than 1.Zhu et al. (2021) previously performed a selective pressure analysis on species from Palaemonidae family within Caridea and also concluded that the Ka/Ks values of all 13 PCGs were less than 1.
The phylogenetic tree results indicated that the genera Plesionika and Heterocarpus within the family Pandalidae were not monophyletic.This finding is consistent with a previous study by Wang et al. (2021) and Liao et al. (2019), who used mitogenome genes or  et al., 2013;Chakraborty et al., 2015;Chakraborty & Kuberan, 2021).This result may be attributed to morphological differences among species within the genus Plesionika, especially the distinct asymmetry of the second pereiopod in P. izumiae and P. lophotes compared to other Plesionika species.At the family level, Pandalidae exhibits strong monophyly, which is consistent with previous research findings (Sun et al., 2020;Ye et al., 2021;Wang et al., 2021;Chak, Barden & Baeza, 2020 Wang, 2021;Cronin, Jones & Antonio, 2022).Regarding phylogenetic relationships among families, the families Acanthephyridae, Oplophoridae, and Alvinocarididae are closely related and form sister groups; similarly, the families Alpheidae and Palaemonidae are also closely related and form sister groups.These results are consistent with previous phylogenetic studies (Chak, Barden & Baeza, 2020;Sun, Sha & Wang, 2021;Cronin, Jones & Antonio, 2022).While our phylogenetic tree topology is consistent with previous research, there are some differences.Specifically, a conflict is noted between our results and those suggested by Li et al. (2011), who proposed that Atyidae represent basal lineages within Caridea, a conclusion based on the analysis of five nuclear genes.Similarly, an inference was made by Bracken, De Grave & Felder (2009) that Atyidae constitute basal lineages within Caridea, an assertion based on the analysis of both mitochondrial and nuclear genes.Furthermore, it was found by Li et al. (2011) that members of the families Palaemonidae and Hippolytidae did not constitute monophyletic groups.The study suggested that members of Hymenoceridae and Gnathophyllidae were clustered within the Palaemonidae clade and that a close relationship was shown between Lysmata amboinensis of Hippolytidae and Janicea antiguensis of Barbouriidae.However, according to the latest records from WoRMS, both Hymenoceridae and Gnathophyllidae have been updated to Palaemonidae (De Grave, Fransen & Page, 2015), and L. amboinensis has also been corrected from Hippolytidae to Lysmatidae (De Grave & Fransen, 2011).Therefore, our study supports that Palaemonidae and Hippolytidae are monophyletic groups.This highlights the importance of incorporating molecular techniques into species identification and classification, as they demonstrate the limitations of past morphology-based species taxonomy.
It has been suggested by some scholars in previous studies that mitochondrial gene arrangement could serve as a new molecular marker to assist in phylogenetic analysis (Zhang et al., 2019;Tan et al., 2018;Wang, Huang & Zou, 2019).Some scholars have also recognized the potential of mitochondrial rearrangement as a ''super'' feature for estimating arthropod phylogenetic (Boore, Lavrov & Brown, 1998;Dowton & Austin, 1999;Dowton, Castro & Austin, 2002;Tan et al., 2017).The relationship between mitochondrial gene arrangement and phylogenetics in Caridea was further analyzed in our study.The phylogenetic tree showed a very clear correlation between both gene arrangement and phylogenetic, and the species having gene rearrangement within each family were clustered together.In the Pandalidae family, the newly sequenced species P. izumiae and P. lophotes exhibited the same gene rearrangement and were closely related in the tree, while three Plesionika species (P.sindoi, P. edwardsii, and P. ortmanni) with the ancestral gene order formed a separate branch.This suggests that the polyphyly phenomenon of Plesionika may be associated with differences in gene order.In Alpheidae, both the Alpheus and Leptalpheus (Williams, 1965) genus underwent gene rearrangement and were clustered together with high support.In the Palaemonidae family, species of the Palaemon genus were clustered together.Except for Palaemon modestus, the same gene rearrangement was undergone by all other species within Palaemon.Conversely, H. picta, experiencing a different gene rearrangement, was positioned alone in a separate cluster.In Lysmatidae, the species with gene rearrangement, L. vittata and E. ensirostris, were also clustered together.Overall, the six families (Pandalidae, Palaemonidae, Alpheidae, Thoridae, Lysmatidae, Hippolytidae) with gene rearrangements were clustered at the base of the Caridea phylogenetic tree, while the five families sharing the same gene order pattern formed a terminal clade.The results of this study show that there is a certain correlation between the phylogenetics of Caridea and the sorting of mitochondrial genes, but additional mitogenomes data are needed to support this result and further investigate their relationship.

CONCLUSIONS
We sequenced the complete mitogenomes of three Plesionika species and analyzed the basic characteristics of these mitogenomes.It was found that these genomes were relatively similar in terms of size, nucleotide composition, and codon usage preference, but exhibited slight structural differences.Additionally, all 13 PCGs in the 12 species of the family Pandalidae underwent purifying selection, with the COI gene experiencing the highest selection pressure, indicating its suitability as an optimal molecular marker for species identification and phylogenetic studies within the Pandalidae.Furthermore, gene rearrangements in the Pandalidae were observed for the first time, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species -P.izumiae and P. lophotes.Phylogenetic analysis revealed a high level of monophyly within the family Pandalidae, but the genus Plesionika appeared to be polyphyletic.By combining the results of gene rearrangements and phylogenetic analysis, a correlation was discovered between the phylogenetics of Caridea and the arrangement of mitochondrial genes.Families that underwent gene rearrangements were located at the base of the Caridea phylogenetic tree, while families without gene rearrangements clustered together at the terminal branch of the phylogenetic tree.This study provides extensive information regarding the mitogenomes of Plesionika, laying a solid foundation for future research on genetic variation, systematic evolution, and breeding of Plesionika using mitogenomes.

Figure 4
Figure 4 The predicted secondary structure of tRNA genes, from trnA to trnW.The nucleotide substitution pattern of tRNA genes in three newly sequenced Plesionika mitogenomes has been exhibited with the reference species P. ortmanni.Full-size DOI: 10.7717/peerj.17314/fig-4

Figure 6
Figure 6 Linear representation of the mitochondrial gene arrangement of the ancestral mitogenome of pancrustaceans and Caridea species.In this study, the three newly sequenced species are marked with blue box.Full-size DOI: 10.7717/peerj.17314/fig-6

Figure 7
Figure 7 The phylogenetic tree based on 13 PCGs was inferred using Bayesian inference (BI) and maximum likelihood (ML) methods.The number at each clade is the bootstrap probability and the three newly sequenced species are marked with red dots.Full-size DOI: 10.7717/peerj.17314/fig-7